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The orbital magnetic moment of a finite piece of matter is expressed in terms of the one-body 
density matrix as a simple trace. We address a macroscopic system, insulating in the bulk, and 
we show that its orbital moment is the sum of a bulk term and a surface term, both extensive. 

The latter only occurs when the transverse conductivity is nonzero and owes to conducting surface 
states. Simulations on a model Hamiltonian validate our theory. 
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According to magnetostatics, the orbital magnetic mo¬ 
ment is defined as m = —dE/d B. This applies to any 
bounded piece of matter. For a homogeneous macroscopic 
system of volume V one writes m = EM, where M is 
the macroscopic magnetization. While in trivial insula¬ 
tors only the bulk states contributes to the magnetic mo¬ 
ment, in nontrivial ones (defined below) a term coming 
from the surface states appears: 

m = m (bulk >+m (surf) , (1) 

and dividing by V we obtain the corresponding contribu¬ 
tions to the orbital magnetization M (bulk) and M (surf) . 

In this paper we propone a new approach to study 
the orbital magnetization, based on an observable which 
allows to discriminate the separate contributions to the 
total magnetic moment coming from the surface and bulk 
states. It applies to any kind of insulator, crystalline or 
noncrystalline. As we will see, a key property of this 
approach is that it is free from the drawbacks related to 
the use of currents. 

We define as “nontrivial” any insulator having a non¬ 
vanishing transverse conductivity at zero B field, which 
we encode in the vector quantity tf 7 = e 7aj 8cr a/ g/2, where 
s 1 af3 is the antisymmetric tensor, and the sum over 
Cartesian indices is implicit. We stress that we are pri¬ 
marily addressing noncrystalline —although macroscopi- 
cally homogeneous—systems; it is nonetheless straight¬ 
forward to extend our treatment to inhomogeneous sys¬ 
tems where M (bulk) varies in space over a macroscopic 
length scale. In any homogeneous insulator the T — 0 
longitudinal conductivity vanishes, while in general a ^ 
0. For a 2D system <x is a dimensionless integer when 
expressed in e 2 /h units, while for a 3D system it has the 
dimension of an inverse length, and is quantized only in 
the crystalline case. 

The difficulties in defining what M really is, and what 
is the role played by the edge states, are closely related 


to the use of currents. Edge states arise due to the ex¬ 
istence of a confining potential but disappear as we con¬ 
sider periodic boundary conditions (PBCs). In the ther¬ 
modynamic limit they do not contribute to the density 
of states per unit volume, but the orbital magnetization 
can be affected by them. However, even in systems where 
only the bulk electrons contribute to the orbital magneti¬ 
zation, the currents which appear at the boundary must 
be taken into account in order to estimate the correct 
value for the magnetic moment. This consideration has 
also been at the root of the longstanding debate about 
the bulk nature of the orbital magnetization.The problem 
is well emphasized in the classic review by Hirst [1]: a fi¬ 
nite magnetized sample is characterized by a dissipation¬ 
less current flowing at its boundary, but in an unbounded 
sample (as addressed in condenser matter physics) the 
macroscopic orbital magnetization M is apparently in¬ 
determinate. In that paper, Hirst analyzes the problem 
in terms of microscopic current densities (either classical 
or quantum), and summarizes the state of the art at the 
time of publication (1997). 

It is clear nowadays that the quantum Hamiltonian 
(and the corresponding ground state) are explicitly 
needed in order to define and to compute M for an un¬ 
bounded sample: the bulk microscopic current density is 
not enough. In fact, it has been shown in 2005-06 [2-4] 
that for a crystalline sample within PBCs: 
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In Eq. (2) BZ is the Brillouin zone, fi is the Fermi level, 
Greek subscripts are Cartesian indices, d a = d/dk a , 
[dk] = dk/(27r) d where d is the dimensionality (either 
2 Or 3), |iijk) = e _ * k r |^k) are the lattice-periodic fac¬ 
tors in the Bloch orbitals, normalized over the unit cell 
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of volume (area in 2D) H ce n; they are eigenfunctions of 
Hk = e~ lk r He lkr , with eigenvalues e^k- 

The existence of a formula for the orbital magnetiza¬ 
tion within PBCs clarifies that, in principle, the orbital 
magnetization can be obtained by considering only the 
bulk of a material. However, the relation between this 
formula and the standard definition of the magnetic mo¬ 
ment within “open” boundary conditions (OBCs), is still 
somewhat obscure. Moreover, the different roles played 
by edge and bulk states are not clearly characterized. 
In fact, the explicit dependence of M from the value of 
H in the bulk, in Eq. (2), implies the presence of edges 
and, possibly, the contact with an external electron reser¬ 
voir in order to change p as a control parameter. It has 
been observed that this is rather puzzling, since Eq. (2) 
addresses a system with no edges [5], whereas any ex¬ 
periment (even gedanken) addresses a bounded sample. 
The approach presented here sheds light on these appar¬ 
ently paradoxical aspects: we provide an alternative ex¬ 
pression for M and, from the very beginning, we get rid 
of currents, which are not the good observable in order 
to analyze the different role played by bulk and surface 
states. In fact, surface currents are due to both bulk and 
surface states and, in general, it is not trivial to disen¬ 
tangle the two contributions. As we will see, in a sense 
a simple analog of our transformation is an integration 
by part: the same integral can be obtained from very 
different integrands. 

The key observation for the following of the present 
work is that Eq. (2) applies as it stands even to a non¬ 
crystalline system, provided we adopt a very large super¬ 
cell, and consequently a mini-BZ. In fact we are adopting 
here the same supercell viewpoint upon which the topo¬ 
logical nature of the quantum Hall effect was established 
[6, 7]. While in any crystalline insulator the spectrum 
e^k is gapped, the spectral gap might close in the large 
supercell limit. This happens for an Anderson insulator, 
where the supercell size must be ideally larger than the 
Anderson localization length. We observe that Eq. (2) re¬ 
tains its validity even for gapless materials [4], and that 
the case where the mini-BZ collapses to a single point 
has been studied in detail [8]. 

We address a macroscopically homogeneous, although 
possibly disordered, piece of matter. For any bounded 
independent-electron system within OBCs, the moment 
is, by definition: 

n N N 
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We recast Eq. (3) in trace form. To this aim we define 
the density matrix (a.k.a. ground state projector) V; we 
will also need the complementary projector Q = I — V . 
Their definitions are 

N oo 

V = J2\ 2= iTb'XTbl- ( 4 ) 

3 =1 j=N+l 

The 7 component of m, from Eq. (3), is 
e ie 

m 1 = -—e ia pTr{Vr a vp} = ~^ £ 7“/5 Tr { Vr a Hrp}, 

(5) 

where we have used v = i[H,r\/h. Lengthy although 
straightforward manipulations of the trace in Eq. (5) 
lead—exploiting antisymmetry— to [9, 10]: 

tc 

m 1 = --^£ ia pTr{Vr a QHQrpV - Qr a VHVr 0 Q}. 

(6) 

We then notice that—when the two traces are expressed 
in the Schrodinger representation—the integrated value 
over the whole sample is the same, but there is a key dif¬ 
ference in the integrands. The unbounded position opera¬ 
tor r is the essential ingredient of Eq. (5), which is there¬ 
fore well defined only within OBCs, i.e. if the | ipj) or¬ 
bitals are square-integrable. Instead, only the projected 
r operator VvQ and its Hermitian conjugate QvV enter 
Eq. (6): this has far reaching consequences. It is known 
since long time that such projected operators are well 
defined and regular even in an unbounded system within 
PBCs. Furthermore (r|'PrQ |r') is nearsighted [11] in in¬ 
sulators, i.e. it decays exponentially (times a polynomial) 
for |r — r'| —> oo. An important theorem proves the ex¬ 
ponential decay even for Anderson insulators [12]. 

Motivated by these considerations, we address the local 
marker in real space for the magnetic moment 

im 7 (r) = ( r\Vr a QHQrpV - Qr a VHVr p Q\r ), 

(?) 

whose integral over the sample gives the total magnetic 
moment m 7 . Thanks to the properties of the operator 
VrQ , the marker ou 7 (r) is well-defined with either OBCs 
or PBCs. Moreover, for an insulator it is local in the bulk, 
since its value in a point r of the bulk is affected only 
by the electronic distribution on a region exponentially 
localized around it. 

We obtain the bulk contribution to the orbital magne¬ 
tization simply by considering the average value of the 
marker (i.e. its integral per unit volume) in the bulk of 
the sample 


We neglect any spin-dependent property here, and 
we deal with “spinless electrons”; v is the quantum- 
mechanical velocity operator, ej, \ifij) are the single¬ 
particle eigenvalues and orbitals, and N is the number 
of electrons. Eq. (3) is the circulation of the whole mi¬ 
croscopic current density: bulk and surface. 


^ bulk) = K> bulk , ( 8 ) 

because the value of im 7 (r) in the bulk is independent of 
the boundary conditions adopted and within PBCs we 
discard any surface effect. Therefore, the surface con¬ 
tribution to the total magnetic moment, m (8urf) =m — 
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m (buik> ^ gj ven by the local marker m 7 (r) on the surface 
of the sample (when OBCs are adopted). As a conse¬ 
quence, the value of m(r) on the surface must be ex¬ 
tensive, although the boundary region is not such. This 
counterintuitive feature has also been confirmed by sim¬ 
ulations at variable sample size (not presented here). 
While the term m ,bulk) is the one ideally measurable by 
accessing the electron distribution in the bulk of the sam¬ 
ple only, m (surf) owes to the electron distribution in the 
boundary region of the sample. As we will see later, 
m ( 3urf) is different from zero only for nontrivial insula¬ 
tors. 

So far, we have implicitly considered an isolated 
bounded sample at fixed N, with the only requirement 
that the resulting Fermi level /i falls in a mobility gap; 
next we are going to consider the same system in contact 
with an electron reservoir which controls the /i (and N) 
value. We observe that the value of M (bulk) (and M (surt) ) 
clearly depends on the (arbitrary) energy zero, while the 
total value of M does not depend on it. It is therefore 
expedient to set the energy zero at the lowest bulk-gap 
edge: with this choice the longitudinal conductivity is 
nonzero for negative /i and vanishes for positive fj, (inso¬ 
far as /i remains in the mobility gap). 

We have defined M (bulk) as a quantity which can be 
computed (and ideally measured) in the bulk of a sample, 
either bounded or unbounded, Eq. ( 8 ); next we wish to 
retrieve M (bulk) within the supercell approach of Eq. (2) 
in reciprocal space. We are going to show that Eq. ( 8 ), 
coincides with Eq. (2), when we set /.t = 0 in the integrand 
only. If localized states in the mobility gap are present 
this still has a /i dependence, as it must be. In order 
to prove the equivalence we need the explicit expression 
for VvQ and its Hermitian conjugate QrV within PBCs. 
While the r operator itself is ill-defined, its off-diagonal 
elements are well defined: this is a staple of linear re¬ 
sponse in solids [13]. One of the known expression is 


all bands, since the difference is a symmetric tensor. Ex¬ 
ploiting completeness we arrive at 


Tr {PrpQHQr a V} 

= Kell y /[dk] (dpUjk\ H k \d a u jk ). (12) 
e Jk <n jBZ 

Similar manipulations performed on the second term 
in the trace in Eq. ( 6 ) lead to the final result 


le 
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By comparing Eqs. (2) and (13) to Eq. (1), we clearly get 
M (surf) by difference: 


M^ urf) = y.^-e-ta/3 "52 /[dk] (<9 a u jk | dpUjy ,). (14) 

hC Sj k<^ BZ 

For cristalline insulators, it can be written as 

= ( is ) 
where C 7 is the Cliern invariant, defined as 

Cry ^ C 7 a/3 ^ ' / [dk] (Oalljk | dpUjk) . (lb) 

ZlT F 

Ej k \ 


However, for a generic insulator, by using Eq. (10) for 
QrV (and its conjugate) we can write Eq. (14) in real 
space [14] as the average in the bulk of another local 
marker c 7 (r) 

< 17 > 

with 


(Vb'kl r \ip jk >) = i^-5(k - k')(uj/ k |V k Wj k ), j ± f, 

(9) 

where we stress that a supercell viewpoint is adopted 
here. We may therefore express QrV as 

QrV = iV, ceii EE/ [d k] iV’i'ic) (ufu |V k Uj k ) {ipjk |. 

ejk<n £j'k>r BZ 

( 10 ) 

The first of the operators entering the trace in Eq. (6) 
becomes thus 

VrpQHQr a V = V C eii E E /["]» 

Ejk<n BZ 

^ d';yk ) (IdpUyk | k ) k (Uy' k d (v: 7/.yk ) (ti/k ■ (H) 

When taking the trace per cell, we may replace the sum 
over the conduction bands (£j' k > /r) with the sum over 


c 7 ( r ) = l~ £ i«0 {r | Vr a QrpV - Qr a VrpQ \ r). (18) 

Notice that Eq. (17), at variance with Eq. (15), is well 
defined both within OBCs and PBCs, and it has the same 
value regardless of the boundary conditions adopted. It 
allows to connect the surface term to transverse conduc¬ 
tivity for a generic insulator. Starting from the standard 
Kubo-Greenwood formula for the transverse conductiv¬ 
ity of a system which in the bulk has a mobility gap, 
straightforward manipulations allow to express even <x 
in terms of VrQ and its Hermitian conjugate QrV [14], 
and we obtain: 

5- = -- M (aurf) . (19) 

[l 

Therefore, as anticipated, M (surf) is different from zero 
for, and only for, nontrivial insulators, that is for insu¬ 
lators having non-zero transverse conductivity. Instead, 
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for trivial insulators M (bulk) = M. While <x is obviously a 
bulk property, M (surf) only depends on something which 
“happens” near its boundary. We stress the virtue of our 
approach: we are not dealing with boundary currents; 
we deal instead with the integrated moment due to the 
boundary states altogether. Eq. (19) is an outstanding 
manifestation of bulk-boundary correspondence: in in¬ 
sulators with nonzero transverse conductivity, the bulk 
and the boundary are “locked”. What appears as “bulk” 
within PBCs, becomes indeed “surface” when addressing 
a bounded sample within OBCs. For this reason the sur¬ 
face contribution to the total M can be “smeared” into 
the bulk of the sample, simply replacing H with H — n 
in the expression for the local marker, Eq. (7): thus even 
the term which is actually due to the surface states of 
the finite sample appears “as if” it were a bulk term. We 
thus recover the local formula for orbital magnetization 
proposed in Ref. [10]. 

Finally, we analyze the variation with // in the gap of 
the quantities M lbulk) and M (aurf) . Here we stress the 
key difference between a system having a spectral gap 
(such as a crystalline insulator) and one having only a 
mobility gap (such as an Anderson insulator). When 
the system is in contact with an electron reservoir, a /i 
variation cannot affect the bulk electron distribution in 
the former case, while the opposite occurs in the latter. 
Therefore, we have a first result: as anticipated, M (bulk) 
is constant with respect to // for crystalline insulators, 
whereas, in the most general case considered here, M (bulk) 
is //-dependent via the projectors V and Q, Eq. (7). 

In order to analyze the variation of M ,surf> with /t, we 
exploit the general relationship: 


For crystalline insulators Eq. (20) reduces to <x = 
—ecd M/9/t, while if there are localized states in the mo¬ 
bility gap the equality no longer holds, and the extra 
term coming from <9M (bulk) /<9/i has to be discounted: as 
shown in Ref. [15] this term corresponds to a current 
non measurable in a transport experiment. By compar¬ 
ing Eq. (20) with Eq. (19), it is immediate to see that 
M< surf > depends linearly on /t. This is a general result, 
valid for any kind of insulator (crystalline or not). 

Simulations addressing model Anderson insulators are 
notoriously very demanding [16]. We can address here 
only a model system having a spectral gap: we therefore 
choose a 2D flake cut from a crystalline system. We adopt 
the Haldane Hamiltonian [17], also adopted by many pa¬ 
pers including Refs. [3, 4, 10, 14], The model is com¬ 
prised of a 2D honeycomb lattice with two tight-binding 
sites per primitive cell with site energy difference of 
2A, real first-neighbor hoppings t\, and complex second- 
neighbor hoppings According to the parameter 

values, the material may have a z = —e 2 /h,0,e 2 /h. All 
of our simulations are performed for rectangular flakes 


with iV Bite =3660 sites, within OBCs. We choose two rep¬ 
resentative cases; their density of states is shown in Fig. 
1 and their bulk magnetization, AH bulk) = TO], bulk) /A, is 
0.037 and 1.006 respectively, in units of tie/(he) (where 
A is the flake area). 

In a tight-binding model the total trace in Eq. (6) be¬ 
comes a discrete sum over the atomic sites of the local 
marker defined in Eq. (7). For the two samples, we show 
in Fig. 2 the local marker value for each site, divided by 
A aite = A/N , ite . The average of this quantity returns the 
orbital magnetization M z . These site contributions are 
not gauge invariant, individually. Only the trace per unit 
area of the local marker is gauge-invariant and in princi¬ 
ple measurable (see below): in this tight-binding model 
it obtains from the sum of any two nearest-neighbor con¬ 
tributions in the bulk region, divided by two. We have 
plotted Fig. 2 for both the trivial insulator with an ar¬ 
bitrary )jl value in the gap, and for the nontrivial one at 
// = 0, i .e. at the bottom of the bulk gap. The absence 
of populated surface states manifests itself in a uniform 
value of the local marker value. The presence of finite 
size effects—due to the discreteness of the spectrum— 
explains the very small spurious boundary contribution 
in the nontrivial case, magnified by the chosen energy 
scale. 

Since we are considering a crystal, in the normal case 
a variation of // in the gap cannot change anything. That 
is why, for the trivial insulator case, Fig. 2 refers to an 
arbitrary // value in the gap. In the nontrivial case, in¬ 
stead, the filling of the boundary states provides a very 
large additional contribution m\ f urf) which scales linearly 
with //,, whereas the value in the bulk remains unchanged. 
Our simulations confirm all these findings. The surface 
nature of the additional contribution is perspicuous in 
Fig. 3 (notice the different scales in Figs. 2 and 3). 

Our theory has addressed magnetization “itself” in 
macroscopically homogenous systems, either disordered 



FIG. 1. (color online) Density of states for the Haldane model 
within OBCs (arbitrary units); the shaded area is the PBCs 
gap. In both cases t 2 = ti/3. Left-hand panel. A = 4t2, 4> = 
0.17r: trivial insulator with a z = 0. Right-hand panel. A = 
t2,4> = 0.4-7/: nontrivial insulator with d z = — e 2 /h. The 
nonzero value in the bulk gap owes to states at the boundary 
of the flake. The vertical dashed line marks the value /.i = 
0.97H (see Fig. 3). 
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C = 0 



FIG. 2. (color online) Plot of im z (r)/ J 4 aite , where im,(r) is 
the local marker for the orbital magnetic moment along z, 
defined in Eq. (7), and A site is the sample area divided by 
the number of sites. The plot is in units of tie/ (he). Top 
panel. Trivial insulator, p in the gap region. Bottom panel. 
Nontrivial insulator, p = 0 at the lowest edge of the bulk gap. 



FIG. 3. (color online) Nontrivial insulator: plot of the differ¬ 
ence Aim z (r)/j4, ite between the cases p = 0.97ti (dashed line 
in Fig. 1) and p = 0. Same units as in Fig. 2. 


or crystalline, while instead—as discussed in Ref. [1] 
one most often measures differences (or derivatives) in 
inhomogeneous situations. Another virtue of the present 
approach is that it applies to inhomogeneous systems 
as well, where M varies on a macroscopic length scale. 
The local marker in the real space is defined in the 
same way, through the diagonal of the relevant operators 
within Schrodinger representation: then the average per 
unit volume becomes the “macroscopic average” of that 
marker, defined as in electrostatics [18]. In the special 
case of an heterojunction the system is insulating when 
p is in the common bulk gap of the two materials. Then 
Eq. (8) yields the bulk magnetizations, while Eq. (17) 
accounts for a /z-dependent interface term in nontrivial 


materials. 

In conclusion, we show that the magnetic moment of 
a macroscopic piece of insulating matter is the sum two 
terms, both extensive, due to states in the bulk and at 
the boundary of the system, respectively, and localized 
in the corresponding regions of the sample. The surface 
term only occurs when the the transverse conductivity is 
nonzero. We stress that we are not addressing the cur¬ 
rent carried by these states: our approach—not based 
on currents—directly yields the moment due to bulk and 
surface states, separately. The approach presented has 
a clear connection with the basic definition of magnetic 
moment for finite systems and analyzes, with a common 
formalism, both crystalline and non crystalline insula¬ 
tors. We have illustrated our theory with a simulation 
for a crystalline system; nonetheless, as said, the theory 
applies as well to an insulator without a spectral gap. 
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